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ABSTRACT 

Observations of gamma ray bursts (GRBs) with Swift produced the initially 
surprising result that many bursts have large X-ray flares superimposed on the 
underlying afterglow. The flares were sometimes intense, had rapid rise and decay 
phases, and occurred late relative to the "prompt" phase. One remarkable flare 
was observed by XRT with a flux >500x the afterglow and a fluence comparable 
to the prompt GRB fluence. Many GRBs have several flares, which are sometimes 
distinct and sometimes overlapping. Short, intense, repetitive, and late flaring 
can be most easily understood within the context of the standard fireball model 
with the internal engine that powers the prompt GRB emission in an active state 
at late times. However, other models for flares have been proposed. The origin 
of the flares can be investigated by comparing the spectra during the flares to 
those of the afterglow and the initial prompt emission. In this work, we have 
analyzed all significant X-ray flares from the first 110 GRBs observed by Swift. 
From this sample 33 GRBs were found to have significant X-ray flares, with 77 
flares that were detected above the 3a level. In addition to temporal analysis 
presented in a companion paper, a variety of spectral models have been fit to 
each flare. A portion of the X-ray flares had enough counts to allow complex 
spectral models, such as Band functions, to be fit to the data. In some cases, 
we find that the spectral fits favor a Band function model, which is more akin to 



department of Astronomy & Astrophysics, Pennsylvania State University, University Park, PA 16802, 
USA 

2 corresponding author email: afalcone@astro.psu.edu 

3 INAF - Osservatorio Astronomico di Brera, Merate, Italy 

4 Universita degli studi di Milano-Bicocca, Dipartimento di Fisica, Milano, Italy 

5 Department of Physics & Astronomy, University of Leicester, Leicester LEI 7RH, UK 

6 NASA/Goddard Space Flight Center, Greenbelt, MD 



-2- 



the prompt emission than to that of an afterglow. While some flares do release 
approximately as much energy as the "prompt" GRB emission, we find that the 
average fluence of the flares is 2.4 x 10~ 7 erg cm -2 in the 0.2-10 keV energy 
band, which is approximately a factor of ten below the average prompt GRB 
emission fluence. We also find that the peak energy of the observed flares is 
typically in the soft X-ray band, as one should expect due to the selection of the 
sample from X-ray data. These results, when combined with those presented in 
the companion paper on temporal properties of flares, supports the hypothesis 
that most X-ray flares are late-time activity of the internal engine that spawned 
the initial GRB; not an afterglow related effect. 

Subject headings: gamma ray bursts, relativistic jets, X-rays, flares 



Introduction 



Since its launch on 2004 November 20, Swift (IGehrels et al. 1 120041 ) has provided de- 
tailed measurements of numerous GRBs and their afterglows with unprecedented reaction 
times. As of 2006 January 24, 110 bursts were detected by the Burst Alert Telescope (BAT; 
Barthelmy et al. 2004). Approximately 93% of these were observed by the narrow field 
instruments in less than 200 ks, and most of those were detected within 200 s (typical reac- 
tion time was ~100 s, but occasionally the BAT detected a burst that was observationally 
constrained). The narrow field instruments are the X-ray telescope (XRT; Burrows et al. 
2005) and the Ultraviolet-Optical Telescope (UVOT; Roming et al. 2005). By detecting 
burst afterglows promptly, and with high sensitivity, the properties of the early afterglow 
and extended prompt emission can be studied in detail for the first time. This also facilitates 
studies of the transition between the prompt emission and the afterglow. The rapid response 
of the pointed X-ray Telescope (XRT) instrument on Swift has led to the discovery that large 
X-ray flares are common in GRBs and occur at times well after the initial prompt emission. 
This paper provides the first survey of the spectral features of a large sample of these X-ray 
flares. 

While there are still many unknown factors related to the mechanisms that produce 
GRB emission, the most commonly accepted mode l is that of a relativisti cally expanding 
fireball with associated internal and external shocks (IMeszaros fc Rees 1119971 ) . In this model, 
internal shocks produce the prompt GRB emission. Observationally, t his emission typicall y 
has a timescale of ~ 20 s for long bursts and ~0.2 s for short bursts (jMeegan et al. Ill996l ). 



The expanding fireball then shocks the ambient material to produce a broadband afterglow 
that decays quickly (typically as with a ~ 1.2 for the nominal afterglow phase). When 
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the Doppler boosting angle of this decelerating fireball exceeds the opening angle of the jet 



into which it is expanding, a steepening of the lightcurve (jet break) is predicted (IRhoads 
19991 ) . For a descript ions of the t heore ti cal models of GRB emissio n and associat e d observa- 



tional prope rties, see [M cszaros | J2002h . Izhang fc Meszaros I (bo04h . Ipiran I (bp05h . Iwbosley 



( 119931 ) , and I Van Paradijs et al. 



the oy erall X-ray lightcurve, s ee 
(bood ). and lwillingalel J2007h . 



For descript i ons of the observatio n al properties of 



Zhang et al. I (l2006ah . iNousek et al. I (l2006luO'Brien et al. 



Several authors have suggested reasons to expect continued ac tivity f r om th e internal en- 
gine of the GRB after the classical "prompt" emission time frame. iKatz I (119971 ) considered a 
model in which a magnetized disk around a central black hole could lead to continued energy 
release in the form of internal shocks. The parameters of this energy release would depend 
on the complex configuration of the magnetic field and the magnetic reco nnection dynamics. 



but ti me periods as long as days for the delayed emission were predicted. iProga and Zhang 



(120051 ) have speculated that energy release can be repeatedly stopped and restar ted at late 
times by magnetic flux accumulation and subsequent release. iPerna et al. I (120061 ) have sug- 
gested that the flares from both short and long bursts can be explained within the context 
of evolution and fragmentation of a viscous accretion disk. For short bursts, in particular, 
Dai et al. I (120061 ) have suggested that late flares can be explained by magnetic reconnection 
events driven by the breakout of magnetic fields from the surface of differentially rotating mil- 



l isecon d pulsars, which resulted from a progenitor compact binary star merger. King et al. 



(120051 ) have speculated that episodic accretion processes could explain continued internal 
engine activity. These authors expect that fragmentation and subsequent accretion during 
the collapse of a rapidly rotating stellar core could explain observations of extended prompt 
emissio n. In general, the domin ant model of an expanding fireball with internal/external 
shocks (IMeszaros fc Rees 1 119971 ) allows for continued prompt emission, provided that the 
internal engine is capable of continuing the energy injection. 

A few observations prior to S wift have included ind ications of flaring from GRBs after 
the prompt GRB emission phase. IWatson et al. I (120031 ) used XMM-Newton to detect line 
emission from GRB 030227 nearly 20 hours after the prompt burst. They inferred continued 
energy injection at this late time, and concluded that a nearly simultaneous supernova and 
GR B event would requ ire sporadic power output with a luminosity in excess of ~ 5 x 10 46 erg 



Piro et al. I (120051 ) used Beppo-SAX to observe two GRBs with relatively small X-ray 



flares. The X-ray flare times for GRB 011121 and GRB 011211 were reported as t=240 s 
and t=600 s, respectively. The spectral parameters of these two X-ray flares were consistent 
with afterglow pa rameters, and these flares were interpreted as the onset of the afterglow 
( Piro et al. I 2005). Two other exa mples of flaring and / or late timescale emission can be 
found in lin't Zand et al. I (120031 ) and iGalli fc Piro I (120061 ) . Although not a detection of late 
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flares from a particular GRB, the work of IConnaughton I (120021 ). in which an ensemble of 
GRBs was analyzed, should also be mentioned. In this study, 400 long GRBs detected by 
the Burst and Transient Source Experiment (BATSE) were analyzed together in the form of 
a summed lightcurve above 20 keV. Significant emission was found at late times (at least to 
1000 s). There are several possible explanations for this emission that do not require flares, 
but flares at various times are certainly one possible explanation. 



Burrows et al. I (j2005bl ) provided the initial report that two bursts detected by Swift 
showed strong X-ray flares. The first of these, XRF 050406, was an X-ray flash with a short, 
and relatively weak, X-ray flare that peaked 213 s after To of the prompt GRB emission. 
Due to the fast rise/decay, the most natural explanation for this flare is continued internal 
engine ac tivity at late times (i. e. delayed prompt emission). XRF 050406 was analyzed in 
detail by iRomano et al. ( 2006 ). Another burst, GRB 050502B, was studied in detail by 
Falcone et al~l J2006h since it was the first dramatic, high-fluence X-ray flare detected. This 
flare, which peaked 740 s after the prompt GRB emission, released as much energy in the 
X-ray band as the prompt GRB released in the 15-150 keV band. Following these two GRBs 
with flares, it became clear that this was a common feature of GRBs, as more and more 
Swift bursts displayed X-ray flares. Although there are a few interesting cases of optical 
flares and of flares simultaneous with higher energy emission detected by Swift-B AT , most of 
these X-ray flares were generally not accompanied by either optical or 15-150 keV emission 
at a detectable level. This implied that the peak of the emission was generally in the soft 
X-ray band. 



There have been severa 


papers studying individual GRBs with X-ray flares ( 


Burrows et al. 


2005b: 


Falcone et a 




2006: 


Romano et a 




2006: 


Paeani et al. 


2006; 


Cusumano et al. 


2006 




Zhane et al. 


2006a 


> 


Vlorris et al. 


2007; 


Goad et al. 


2007; 


Krimm et al. 1 


2007 


). While the 



detailed study of individual flares is important, it is equally important to look at the proper- 
ties of the flares in a more general sense to look for general trends and overall mean properties 
of the flares. By comparing these overall properties to those of the prompt GRB emission 
and the afterglow emission, the mechanism of the flare emission may be elucidated. Further- 
more, we can see if there are multiple classes of flares, or if the flare parameters all fall into 
one uniform distribution. 



In this paper and a companion paper (jChincarini et al. 1 120071 ) (referred to as Paper I 
hereafter), we present the first temporal and spectral study of a statistical sample of X-ray 
flares within GRBs. Paper I presents the temporal properties of the sample, and this paper 
presents the spectral properties of the sample. The sample includes all bursts, up until 2006 
Jan 24, for which Swift detected at least one significant X-ray flare. 
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2. The Sample 

The initial sample was chosen by looking at all Swift-XBT light curves, between launch 
and 2006 January 24, and eliminating the ones that did not show any hint of deviation from 



a power law de c ay wit h typical breaks (see IZhang et al. I (j2006al ) , iNousek et al. I (120061 ) , and 



O'Brien et al. I (120061 ) for a discussion of typical lightcurve breaks). The remaining light 
curves were then fit with a broken power law decay (referred to as the underlying decay 
curve and subscripted throughout text with 11 UL") superposed with a power law rise and 
decay for any flares that appeared above this underlying decay curve. The start time of 
the flare was then defined as the time that the power law rise of the flare intersected the 
underlying decay power law. Similarly, the stop time of the flare was then defined as the 
time that the power law decay of the flare intersected the underlying decay power law. These 
times, t s tart and t stop , are defined relative to the trigger time, T , of the GRB. The signal to 
noise of the flare was then defined, using simple Poisson statistics, as: 

S/N= (1) 

VNtotal+NuL 

where N tota i is defined as the total number of photons during the flare time interval, and Njjl 
is defined as the number of photons from the fit to the underlying decay curve during the 
flare time interval. Only the flares with S/N> 3 were retained in the sample. This analysis 
of 110 GRBs resulted in 33 GRBs with at least one significant flare, and it resulted in a total 
of 77 flare time intervals, which are listed in Table [H Some of these 77 time intervals overlap 
one another so it is not always clear where one flare begins and another ends. We define 
the start and end of each flare as described in the analysis section below. This sample and 
the sample of Paper I are largely overlapping, but they differ somewhat due to the different 
approach and goals. The flares for which temporal properties can be obtained are frequently 
different from those for which spectral properties can be obtained. 
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Table 1. The Flare Sample 



GRB 


Flare 


t-start (s) 


tstop (s) 




S/N 


GRB050219a 


1 


118 


453 


120 


18.5 


GRB050406 


1 


139 


361 


205 


11.3 


GRB050421 


1 


136 


165 


156 


3.4 


GRB050502b 


1 


410 


1045 


695 


145.7 


GRB050502b 


2 


19958 


48591 


29896 


7.2 


GRB050502b 


3 


50457 


178280 


75355 


18.4 


GRB050607 


1 


94 


255 


145 


10.3 


GRB050607 


2 


255 


640 


312 


25.2 


GRB050712 


1 


88 


564 


252 


31.0 


GRB050712 


2 


302 


435 


339 


12.9 


GRB050712 


3 


415 


590 


478 


8.9 


GRB050712 


4 


788 


952 


888 


3.8 


GRB050713a 


1 


101 


155 





11.7 


GRB050713a 


2 


155 


210 





3.2 


GRB050714b 


1 


285 


832 


374 


19.2 


GRB050716 


1 


155 


211 


177 


11.2 


GRB050716 


2 


315 


483 


385 


13.2 


GRB050724 


1 


78 


230 


120 


102.6 


GRB050724 


2 


63 


342 


261 


33.7 


GRB050724 


3 


13406 


402320 


55783 


19.7 


GRB050726 


1 


151 


195 


162 


3.0 


GRB050726 


2 


219 


324 


274 


12.2 


GRB050730 


1 


210 


280 


228 


20.6 


GRB050730 


2 


323 


611 


435 


51.9 


GRB050730 


3 


611 


795 


678 


33.7 


GRB050730 


4 


9654 


12578 


10319 


33.2 


GRB050802 


1 


312 


457 


435 


3.8 


GRB050803 


1 


513 


879 


753 


5.8 


GRB050803 


2 


889 


1516 


1116 


4.3 


GRB050803 


3 


4455 


5703 


5367 


5.8 


GRB050803 


4 


7345 


27698 


22669 


14.2 


GRB050803 


5 


7646 


13093 


11613 


14.0 


GRB050803 


6 


17240 


27698 


18873 


5.1 


GRB050814 


1 


1133 


1974 


1350 


3.0 


GRB050814 


2 


1633 


2577 


2138 


6.1 


GRB050819 


1 


56 


253 


174 


11.5 


GRB050819 


2 


9094 


36722 


19733 


6.2 


GRB050820a 


1 


200 


382 


234 


66.6 


GRB050822 


1 


106 


190 


143 


21.3 


GRB050822 


2 


212 


276 


240 


8.4 


GRB050822 


3 


390 


758 


433 


50.9 


GRB050904 


1 


343 


570 


463 


41.6 


GRB050904 


2 


857 


1141 


953 


3.0 


GRB050904 


3 


1149 


1343 


1235 


4.2 


GRB050904 


4 


5085 


9001 


6765 


23.0 


GRB050904 


5 


16153 


24866 


17329 


22.1 
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3. Analysis 

The data were reduced using the latest HEAsoft tools (version 6.1.0), including Swift 
software version 2.0, and the latest response (version 8) and ancillary response files (created 
using xrtmkarf) available in CALDB at the time of analysis. Data were screened with 
standard parameters, including the elimination of time periods when the CCD temperature 
was warmer than -48° C. When analyzing WT data, only grades 0-2 were included, and 
when using PC mode data, only grades 0-12 were included. Source and background regions 
were both chosen in a way that avoids overlap with serendipitous sources in the image. For 
PC mode data, the source region was typically a 20 pixel (47 arcsec) radius circle. The 
background region was typically a circle with a radius of 60 pixels chosen in a source-free 
region (40 pixels if the field was crowded). All quoted errors are la unless otherwise stated. 

In order to avoid pile-up effects in some of the higher count rate PC mode data (> 
0.5 c/s), an annular source extraction region was used with an inner radius that varied as a 
function of rate. WT mode data is free of significant pile-up effects for nearly all of the flares. 
However, pile-up does begin to have a marginal systematic effect on WT mode data above 
100 count/s. For the few flares in this sample with a brief excursion above 100 count/s, the 
effect of pile-up is insignificant on this analysis since it averages the spectrum over the entire 
time interval of the flare, and as a result, the vast majority of flare photons are not affected 
by pile-up. 

3.1. Light Curve Analysis 

The light curves, and the corresponding temporal analysis, are presented in Paper I. 
However, it is worth mentioning a few of the temporal analysis issues related to the spectral 
analysis presented here. In particular, there are a few differences between the approaches 
of this paper and Paper I. The t start and t stop times in our analysis were not constrained 
to be exactly the same as those of Paper I. However, they are approximately the same and 
the differences are irrelevant for the purposes of fitting the spectra. The small differences 
arise from the fact that this analysis fits temporal power law curves to the rise and decay 
portions of the flares, whereas Paper I fits Gaussians to the flares. The points on the light 
curve where these power laws intersect the underlying decay curve power law are defined 
as t s tart and t st0 p, and they are reported in Table [H This method allows us to easily define 
a temporal region for performing spectral fits on flare data, even if the flare is missing a 
large fraction of its lightcurve for any reason. We are still able to fit a spectrum to large 
flares that are missing some data on the rising or falling portion of the flare, even if they 
do not have a well constrained temporal fit. When calculating the fluence in that portion 
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Table 1 — Continued 



GRB 


Flare 


i^start \p) 


t t M 

L stop \ & } 


t , M 

''peak V J 


S /N 


GRB050904 


6 


18383 


38613 


24156 


19.5 


GRB050904 


7 


25618 


30978 


29392 


21.6 


GRB050908 


1 


129 


306 


145 


7.3 


GRB050908 


2 


339 


944 


404 


14.0 


GRB050915a 


1 


55 


170 


111 


14.3 


GRB050916 


1 


16755 


32357 


18898 


20.1 


GRB050922b 


1 


357 


435 


377 


12.1 


GRB050922b 


2 


476 


560 


497 


5.6 


GRB050922b 


3 


630 


1541 


827 


39.4 


GRB051006 


1 


115 


148 


132 


9.6 


GRB051006 


2 


132 


201 


162 


7.5 


GRB051006 


3 


330 


749 


495 


7.5 


GRB051016b 


1 


374 


1940 


483 


3.1 


GRB051117a 


1 


2 


4322 


157 


117.6 


GRB051117a 


2 


134 


2794 


380 


124.1 


GRB051117a 


3 


292 


1313 


628 


70.4 


GRB051117a 


4 


574 


2695 


926 


78.6 


GRB051117a 


5 


642 


1820 


1097 


71.0 


GRB051117a 


6 


1237 


3119 


1335 


95.3 


GRB051117a 


7 


659 


3126 


1535 


85.8 


GRB051210 


1 


115 


152 


132 


4.4 


GRB051227 


1 


86 


245 


120 


16.7 


GRB060108 


1 


193 


429 


285 


2.1 


GRB060108 


2 


4951 


37986 


10471 


6.1 


GRB060109 


1 


4305 


6740 


4810 


5.0 


GRB060111a 


1 


27 


196 


110 


51.5 


GRB060111a 


2 


109 


203 


171 


38.6 


GRB060111a 


3 


215 


433 


312 


107.4 


GRB060115 


1 


331 


680 


406 


8.7 


GRB060124 


1 


283 


644 


574 


222.6 


GRB060124 


2 


644 


1007 


694 


179.8 
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of the flare, a correction factor is applied based on a power law extrapolation of the flare 
light curve. To define the underlying decay curve, we use multiply broken p ower laws that 



account for the various phases of the GRB and afterglow light curve decays (INousek et al. 



20061 : IZhang et al. Il2006al : lO'Brien et al. 1120061 ) 



In some cases, the time range for spectral extraction did not include the entire flare 
time range due to reasons such as incomplete light curves or overlapping flares. The time 
regions used for spectral extraction are shown in Table [2J This table shows the times for 
flare spectral extraction and underlying light curve spectral extraction. In some cases, the 
underlying spectral extraction used multiple time regions to improve statistics. In other 
cases, there was one large contiguous time period for underlying spectral extraction. 
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Table 2. The time regions used for spectral extraction of flares and the time regions used 
for the extraction of underlying light curve spectra for each flare are tabulated. In many 

cases, the underlying spectra was constrained with one time region with sufficient photons 
to obtain spectral parameters, and thus there are no entries in the last four columns. In 
some cases, statistics were maximized by using multiple time regions for the underlying 

portion. In a few instances there are dashes for all underlying time regions, indicating that 
the canonical value for the underlying spectral index was used, as described in the text. 



GRB 


Flare 


Flares 
tbeqin t s ) 


tend V s ) 


Underlying 
( 1 )t b (sj 


tend i s ) 


tbegin ( s ) 


tend ( s ) 


tbegin t s ) 


tend ( s ) 


GRB050219a 


1 


118 


453 


670 


29603 












X 


139 


361 


1447 


ij ±iJOO\J 










GRB050421 




136 


165 


167 


488 










GRB050502b 


x 


410 


1045 


5384 


20369 


161890 


299820 






GRB050502b 


2 


19958 


48591 


57 


355 


1545 


19958 


178280 


264880 


GR.B050502b 


3 


50457 


178280 


57 


355 


1545 


19958 


178280 


264880 


GRB050607 


x 


94 


255 


92 


94 


685 


20997 






GRB050607 


2 


255 


640 


92 


94 


685 


20997 






GRB050712 


x 


88 


299 


5157 


105060 










GRB050712 


2 


302 


435 


5151 


77682 










GRB050712 


3 


415 


590 


5074 


63858 










GRB050712 


4 


788 


952 


5074 


63858 










GRB050713a 


x 


101 


155 


3541 


399630 










GRB050713a 


2 


155 


210 


3541 


399630 










GRB050714b 


x 


285 


542 


3639 


139690 










GRB050716 


x 


155 


211 


105 


155 


211 


331 






GRB050716 


2 


315 


483 


211 


331 










GRB050724 


x 


78 


230 


433 


27350 










GRB050724 


2 


222 


342 


433 


27350 










GRB050724 


3 


13406 


402320 


433 


27350 










GRB050726 


x 


151 


195 


324 


12646 










GRB050726 


2 


219 


324 


324 


8358 












x 


210 


280 


132 


210 


280 


313 








2 


323 


611 














GRB050730 


3 


611 


795 














GRB050730 


4 


9654 


12578 


4366 


6863 


26422 


99149 






GRB050802 


x 


312 


457 


494 


2873 










GRB050803 


x 


513 


879 


34808 


778510 










GRB050803 


2 


889 


1516 


34808 


778510 










GRB050803 


3 


4455 


5703 


34808 


778510 










GRB050803 


4 


7345 


27698 


34808 


778510 










GRB 05 08 03 




10396 


13093 


34808 


778510 










GRB050803 


Q 


17240 


27698 


34808 


778510 










GRB050814 


x 


1133 


1974 


5646 


8644 


32429 


98328 






GRB050814 


2 


1633 


2577 


5774 


8741 


32794 


96149 






GRB050819 


x 


154 


193 














GRB050819 


2 


9094 


36722 


475 


7975 


36722 


55757 






GRB050820a 


1 


200 


258 


4811 


5099900 










GRB050822 


1 


106 


190 


5692 


4932900 










GRB050822 


2 


212 


276 


5911 


4795400 










GRB050822 


3 


415 


616 


4714 


5628400 










GRB050904 


1 


343 


570 


586 


868 










GRB050904 


2 


857 


1141 


588 


876 










GRB050904 


3 


1149 


1343 


588 


861 










GRB050904 


4 


5085 


7110 


581 


865 










GRB050904 




16153 


18205 


586 


873 










GRB050904 


6 


22221 


25379 


586 


873 










GRB050904 


7 


27854 


30978 


586 


873 










GRB050908 


1 


129 


306 














GRB050908 


2 


339 


944 














GRB050915a 


1 


55 


170 


170 


7424 










GRB050916 


1 


16755 


32357 


221 


13085 










GRB050922b 


1 


357 


435 


348 


355 


435 


476 


560 


623 


GRB050922b 


2 


476 


560 


348 


355 


435 


476 


560 


623 


GRB050922b 


3 


630 


1541 


348 


355 


435 


476 


560 


623 


GRB051006 


1 


115 


148 














GRB051006 


2 


148 


180 














GRB051006 


3 


330 


749 














GRB051016 


1 


374 


194(3 


3778 


382750 










GRB051117a 


1 


113 


231 


16046 


2410600 










GRB051117a 


2 


295 


571 


16046 


2410600 










GRB051117a 


3 


571 


729 


16046 


2410600 










GRB051117a 


4 


817 


1044 


16046 


2410600 










GRB051117a 


5 


1044 


1237 


16046 


2410600 










GRB051117a 


6 


1237 


1466 


16046 


2410600 










GRB051117a 


7 


1466 


1737 


16046 


2410600 
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3.2. Spectral Analysis 



Spectral models were fit to data using Xspec version 12.2.0 (lArnaud et al. II2005I ). Spec- 
tra were fit in the 0.3 to 10.0 keV energy range. A systematic error of 3% was assigned 
throughout the energy range due to uncertainties in the response of the instrument, par- 
ticularly below 0.6 keV. During fitting, all spectra were binned to >20 photon/bin, and \ 2 
statistics were used. 

This work is attempting to apply four different models to flare spectral data. However, 
it is clear that the flare itself is not the only X-ray emission during the time of the flare. The 
underlying afterglow of the GRB is usually already in progress at the time of the flare. In 
some cases, this is a small fraction of the flux from the flare and would merely add a small 
systematic effect to the spectral fit. However, in other cases, the underlying light curve is 
a large fraction of the observed X-ray emission, and its effect must be taken into account. 
We choose to do this by selecting a region of the light curve before and/or after the flare 
and fitting the spectra in this time region to a simple absorbed power law, which has the 
following form. 

fuL = C UL [e- N ^ E \^]- T ^ (2) 
where Nh,ul is the neutral Hydrogen column density wi th units atoms / 'cm~ 2 , a(E) is th e 



energy dependent photoelectric absorption cross section ( Morrison and McCammon II1983I ). 
T is the spectral photon index, and C is the normalization constant in units of photons 
cm -2 s -1 keV" 1 . If the fit results in a value for the Nhjjl that is significantly below the 
Gala ctic Nh, then the fit is rec alculated with Nh,ul set equal to the Galactic value taken 



from iDickey & Lockman I (Il990l ). In cases for which there are not enough photons that are 
obviously part of the underlying light curve (i.e. independent of the flare), a canonical value 
of T = 2.0 is chosen, and the N H is simply tied (i.e. forced to be equal) to the Hydrogen 
column density in the subsequent flare fitting, Nujiare- 



Table 2 — Continued 



GRB 


Flare 


Flares 




Underlying 
(1) * 6e9i7l (s) 


(1) te„ d (S) 


< 2) w» W 


Vt end (s) 




<3) t c „ d (s) 


GRB051210 


1 


115 


152 


162 


426 










GRB051227 


1 


86 


245 


258 


20156 










GRB060108 


1 


193 


429 














GRB060108 


2 


4951 


37986 














GRB060109 


1 


4305 


6740 


8784 


325220 










GRBOeDUla 


1 


75 


137 


2905 


712320 










GRB060111a 


2 


145 


204 


2905 


712320 










GRB060111a 


3 


215 


433 


2905 


712320 










GRB060115 


1 


331 


680 


117 


257 










GRB060124 


1 


283 


644 


10605 


14232 


32067 


74305 






GRB060124 


2 


644 


1007 


10458 


14432 


33443 


71248 
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The normalization of this spectral power law is then found by using the power law 
fit to the temporal decay of the underlying light curve before and/or after the flare. The 
underlying temporal power law is extrapolated into the time region of the flare, and it is 
integrated over that time range to obtain the expected counts from the underlying decay 
during the flare. This provides the scale factor by which Cjjl must be normalized. 

Once the underlying spectra and flux have been estimated, as described above, these 
values are frozen and used as an additive part to the four spectral models for the flares. We 
then attempt to fit this additive model to the overall flare plus underlying afterglow spectral 
data. The four models we apply are: simple power law, exponentially cut off power law, power 
law plus blackbody, and Band function. The application of these non-power-law models, with 
more complex curvature, has been motiva t ed by t he sim il ar ap p lication of models to p rompt 



GRB emission surveys (e.g. Band et al. I (I19931 ): iRvde I (120041 ); Kaneko et al. I (12009 1). For 



all four models, we also apply photoelectric absorption, which is free to vary. To illustrate 
the method, the equations of two of the models are shown below. 



Simple absorbed Power Law: 



total 



c 



flare [ 



L keV 



keV 1 



(3) 



Absorbed Exponentially Cut off Power Law: 



ft 



total 



c 



flare [ 



rp rp 



■keV 1 



keV 1 



(4) 



The other two models used are the thermal blackbody plus power law model and the 
Band function model, both of which are added to the underlying afterglow power law model 
in the same way as the exponential cutoff power law model is added to the underlying model 
i n the equations sh own above. The blackbody model and the GRB Band function model 
( IBand et al. Ill993l ) are described by lArnaud et al. I (120051 ). In all cases, Cul, N H)U l, and 
Ttjl are frozen to the values determined using data from a region before and/or after the 
flare, and all flare parameters (e.g. C f i are , N HJlare , T flare , E 0jflare , E c j lare , kT flare ) are free 
to vary during the fitting process. 



4. Spectral Results 

The results from applying the spectral models described above are presented in this 
section. In order to maximize the prospects for having reasonably constrained parameters, 
we selected only the flares for which there were more than 15 degrees of freedom during the 
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fitting of the power law flare model. From this point forward in the text, we refer to these 
flares as Gold flares. 



4.1. Overall Spectral Parameters of Underlying Decay 

The spectral parameters derived by fitting an absorbed power law model to the under- 
lying afterglow of GRBs with Gold flares are shown in Figure [TJ The mean of the photon 
index distribution is 1.9 with a standard deviation of 0.3. This is consistent with the typical 
photon index for GRB afterglows. 

4.2. Overall Spectral Parameters of Flares 

The spectral parameters derived by fitting the absorbed power law model to the Gold 
flares are shown in Figure [2] and Table [31 The index of each flare, shown in column 1 of each 
table, corresponds to the x-axis of the plots. A simple absorbed power law can provide a 
reasonable fit in most cases. The spectral parameters derived by fitting the absorbed Band 
function model to the Gold flares are shown in Figure [3] and Table HI Once again, the fit 
is reasonable in nearly all cases. The spectral parameters derived by fitting the absorbed 
exponential cutoff model to the Gold flares are shown in Figure H] and Table El and the 
spectral parameters derived by fitting the absorbed power law plus blackbody model to the 
Gold flares are shown in Figure El and Table El 

It is clear that there are many cases for which a power law provides a satisfactory 
fit. However, it is also clear that there are many cases for which a more complex model, 
such as a Band function, provides a superior fit. In order to explore the distribution of the 
change in the quality of fit, a histogram of the A% 2 between the power law fits and the Band 
function fits is shown in Figure El This histogram shows the x\ ow ~ XBand f° r the 47 Gold 
flares. The mean degrees of freedom were 130 and 128 respectively. For comparison, we have 
also simulated this same Ax 2 for a fake distribution of power law spectra. To do this, we 
simulated 1000 fake spectra that were power laws with spectral photon indices 1.9±0.3 and 
the same mean degrees of freedom as the spectra in the flare sample. These spectra were then 
fit in the same way that we fit the spectra of the flare sample, and the A% 2 was calculated. 
The resulting Ax 2 distribution from this simulation was plotted as a curve overlaying the 
histogram of the real data in Figure El 

As expected, there are many flares that can be fit equally well by both models, but some 
flares have a better fit using a complex model such as the Band function. The distribution 
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Fig. 1. — Properties of absorbed power law spectral fits to data from a region of the lightcurve 
in which no flares were present, for all GRBs with Gold flares (i.e. these are the spectral 
parameters of the underlying light curve). The index number of the flares shown on the 
x-axis simply refers to the index number for each flare shown in column 1 of Table [3j 
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Power Law Model Parameters for Flares 
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Fig. 2. — Properties of power law spectral fits to flare data for all Gold flares. The index 
number of the flares shown on the x-axis simply refers to the index number for each flare 
shown in column 1 of Table [HI The top panel corresponds to the fit for the neutral Hydrogen 
column density, the second panel corresponds to the photon index (Tfi are ), and the bottom 
panel is the reduced x 2 f° r eacn fit- 
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Table 3. Properties of power law spectral fits to Gold flares 



Index Number 


GRB 


Flare 


JV H (10 20 cm- 2 ) 


Photon Index 


Xred 


DOF 



10 

11 

12 
13 
14 
15 
16 
17 
18 
19 
20 
21 
22 
23 
24 
25 
26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 



GRB050219 
GRB050502 
GRB050502 
GRB050607 
GRB050712 
GRB050712 
GRB050713 
GRB050713 
GRB050716 
GRB050716 
GRB050724 
GRB050724 
GRB050724 
GRB050726 
GRB050730 
GRB050730 
GRB050730 
GRB050730 
GRB050802 
GRB050803 
GRB050803 
GRB050820 
GRB050822 
GRB050822 
GRB050822 
GRB050904 
GRB050904 
GRB050904 
GRB050904 
GRB050904 
GRB050916 
GRB050922 
GRB050922 
GRB050922 
GRB051117 
GRB051117 
GRB051117 
GRB051117 
GRB051117 
GRB051117 
GRB051117 
GRB051227 
GRB060111 
GRB060111 
GRB060111 
GRB060124 
GRB060124 



39.9 
11.7 
10.7 



+11.1 
9.1 
-0.7 
-0.7 
+ 5.3 
4.5 



22.5 



+ 5.3 



-4.7 
(-4.7 
J -4.3 
20.5+ 89 



24. cr 



52. 1 1 



-6.9 
,+5.9 
-5.3 
50. C_ 10 9 

24 4+ 86 - 4 
z *'*-24.4 

B2 .it?;S 

55.6tt 6 3 
27.4+J 8 ; 8 
14.6t 9 3 3 9 

»-8±i.5 

8 4+ 1 ' 1 
5.9 ' ' 



+ 1.3 
1.2 
+ 3.4 
3.0 
+ 8.8 
6.3 
+ 9.0 
7.8 
+ 28.3 
21.7 
+ 1.6 
-1.5 
+ 7.6 
-6.1 
+4.9 
4.3 
+ 22.0 
-14.4 

14 6 +2 ' 6 
lq - D -2.4 

8 9+ 2 ' 4 
°- y — 2.2 

1 ^- u -3.0 
4 2+ 2 ' 5 
+ 2.5 
2.2 



13.0 
7.7 
47.9 
87.8 
9.5 
11.1 
19.4 
29.9 



7.5"' 



99.4 
47.9 
20.3 
14.6 
16.7 
16.8 
13.5 
14.4 
14.1 



+ 32.6 
26.9 
+ 16.4 
12.6 



+ 18.3 



-12.0 
+ 2.4 
-2.2 
+ 1.1 
1.1 
1.0 
-0.9 
b + 1.5 
-1.4 
+ 1.4 
1.3 
+ 1.5 
1.5 

16.0+ 12 



-1.1 
K+l-2 
-1.2 
a + 9.1 
-7.1 
-+3.9 
-3.6 
,+4.1 
-3.7 
26.5+ 1 ' 4 



12.57 
29. 9j 
38. 5+ 
31. 1" 1 



2 67+ 41 
^°'-0.34 
, ,, + 0.04 
a.oo_ 04 

2 10+ ' 27 
z ' lu -0.23 

2 40+ 023 
z,w -0.20 

2 l 3 +° 18 
z ' 1,3 -0.17 

, ns+ SS9 

■ s - vs -0.44 
2 19+ ' 11 

s-sotgli 

1 22+ '' 93 
- L " s -0.75 
, ,0+0.93 

°-°°-0.63 
, __+0.03 
L - 1 1 -0.03 

2 94+ 13 
z ' a ^-0.12 

1.61 



2 

1.71 



+ 0.15 
0.13 
+ 0.68 
-0.50 



,+0.12 
-0.12 
1 R«+ - 05 
1 ' 66 -0.05 

1 92+ 007 
1 ' az -0.07 

2 20+ ' 14 
z - zu -0.13 
, ,, + 0.46 
^'"-0.36 
2.27+ 



4.55 
0.82 



-1.4 
^+0.5 
-0.5 

16 4 4 " - 5 
1D - -0.5 



18. 4 1 



+ 0.23 
0.21 
+ 1.12 
0.82 
+ 0.04 
0.04 
, 7 „ + 0.30 

2.86±«; 2 I 

, ,,+0.09 

1 96+ 010 
laD -o.io 

1 96+ 014 
1 - aD -0.13 

1 - S1 -0.12 

is& _0.11 
+ 0.31 
0.30 
+ 0.78 
0.58 
+ 0.86 
-0.60 

2 36+ 010 

1 uo + 0.04 
1 ' 88 -0.04 
„ ,,+0.04 
^■^-0.04. 

2 26+ 07 
z - ZD -0.07 

+ 0.06 
0.06 
+ 0.08 
0.08 
+ 0.05 
0.05 
+ 0.06 
0.06 
+ 0.15 
0.14 
+ 0.14 
0.13 

2 ' 86 -0.17 
, , 7 +0.05 

, ,,+0.01 
^ -0.01 

1 67+ ' 02 
1D '-0.02 



1.70 
3.94 
2.66 



2.133 
2.51 4 
2.22+ 
2.25j 
1.53+ 
2.89 4 
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Fig. 3. — Properties of Band function spectral fits to flare data for all Gold flares. The 
index number of the flares shown on the x-axis simply refers to the index number for each 
flare shown in column 1 of Table HI The top panel corresponds to the fit for the neutral 
Hydrogen column density, the second panel corresponds to the low energy photon index (a), 
and the third panel corresponds to the high energy photon index (/?). The fourth panel 
is the e- folding energy (E ), which is related to the peak spectral energy by the relation 
Epeak = (2 + q)Eq. The bottom panel is the reduced \ 2 f° r the fits. 
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Table 4. Properties of Band function spectral fits to Gold flares 



Index Number 


GRB 


Flare 


JV H (10 20 cm- 2 ) 




3 


E peak (keV) 


"red 


DOF 


1 


GRB050219 


1 


30.8+f, 3 

— 3.3 


1.15+?- 43 

— 1.07 


2 -5 2 io27 


0.6+0- 5 . 

— 0.5 


1.12 


36 


2 


GRB050502 


1 


3.3+li 


0.74+°- 11 


2 ,, + 0.'04 
— 0.05 


i.o+;-j 


1.20 


326 


3 


GRB050502 


3 


9.2+fl 


1.06+?-'^ 


2 -0 1+ S'26 


l.O+oooo 


0.85 


29 


4 


GRB050607 


2 


16.7+f/ 


1.11+0-" 


2.31+0-1 8 


J g + 66.8 


0.87 


31 


5 


GRB050712 


1 


28.2+8-0 


1.16+f-ll 


2.08+0-JO 




1.47 


5 5 


6 


GRB050712 


2 


17.2+1 1 ' 5 3 


1.70+011 


2.60+»- 3 « 


0.8+0- 3 


0.98 


16 


7 


GRB050713 


1 


44.1+1° 7 8 


1.69+" 3 ? 

— 0.57 


S.97+8- 87 

— 1.03 


6.5+ 21 „ 9 


1.02 


186 


8 


GRB050713 


2 


36. 1+ 44 , 7 


1.44+«|1 


3.02+0- 34 


0.9+ 42 / 


1.24 


46 


9 


GRB050716 


1 


S6.7+ 97 - 4 
— 3b. 7 


1.02+ 2 - 07 

— 8.74 


1.28+0™ 


2.6+1% 


0.28 


54 


10 


GRB050716 


2 


7.9+11 

— 7.4 


1.23+0^ 


„ -+19.37 
— 9.37 


i.o+i- 9 

2.7+ 2 -'l 

— 1.0 


0.34 


53 


11 


GRB050724 


1 


36.5 +r ; \ 

— 3. 3 


0.69+H 3 


1.87+000 

— 0.07 


1.01 


328 


12 


GRB050724 


2 


43. 6+|| 


1.85+ ,' 99 


3.11+0- 27 


2-3+o1 


0.97 


52 


13 


GRB050724 


3 


13.4+1 3 ; 


0.51+ o: « 2 


1.91+0- 4 ; 


2.7+ 2 - 7 


1.30 


20 


14 


GRB050726 


2 




— .88 


2.92+. O,- 77 


1.2+fof 


0.92 


35 


15 


GRB050730 


1 


3.2+3?, 




1.77+0- 16 


1.8+1' 3 


1.03 


57 


16 


GRB050730 


2 


6.7+\i 


0.63+H 3 

— 0.05 


i.6i+o ; | 

— 0.03 


0.9+01 

— 0.9 


0.91 


185 


17 


GRB050730 


3 


5.2+\l 


0S6+_°ll 


1.80+00 7 


0.9+0- 2 


0.84 


104 


IS 


GRB050730 


4 


0.9+ii 


0.38~^ri r\£r 

— 0.9o 


2.40+0-' 24 

— 0.32 


i.i+;'t 


0.96 


79 


19 


GRB050802 


1 


1.8+]'i 

— 1.8 


1 17 +0.45 
— 0.52 


9.16+1 9 ,! 6 

— 9.16 


2.4+S-l 
— 0. o 


0.97 


29 


20 


GRB050803 


5 


43 - 1 + 22 1 4 


0.97+0-g 


9 9 o + 0.18 
19 


0-2 + o'2 


1.45 


32 


21 


GRB050803 


6 


64.8+ 2 1 4: « 


2.26+0f= 


9.37+1O3 37 


0.3+1 = 


1.17 


16 


22 


GRB050820 


1 


4.7+ 4 - 7 




0.82+0 04 


3.0+ 2 -'! 


1.12 


201 


23 


GRB050822 


1 


4.0+H 


0.90+1-1 4 


1.91+ - 40 


2.6+ 2 -l 

— 2.5 


0.42 


25 


24 


GRB050822 


2 


S.9+1% 1 


1.14+° 13 


2.82+o|| 


0-4+ i 


1.12 


29 


25 


GRB050822 


3 






8.61+1| 9 


0.9+11/ 


1.93 


16 


26 


GRB050904 


1 


13 - 4+ 10 9 2 


0.95~^~q 


, 7r +0.09 
09 




0.98 


180 


27 


GRB050904 


4 


3.3+li 


1.30+0'f, 

— 0.35 


9.18+ 7 '°0 

— 0.18 


4.2+ 2 '- 2 


1.30 


36 


28 


GRB050904 


5 


1 7+ 6 -' 1 
'-1.4 


0.29+° ; 7 ? 


1 89+ 013 


9+ 1:3 


0.89 


24 


29 


GRB050904 


6 


5 0+ 5 - 


1 60+ - 78 
1-D — 0.36 


1 81+013 




1.04 


21 


30 


GRB050904 


7 


4 0+ 4 - 5 




9 , 9 + 0. 38 
zl ^-0.18 


5.8±S;« 


0.93 


22 


31 


GRB050916 


1 


q7 „+60.1 


1 00 + 2 ' 53 


, K o + 0.22 
l' bS -0.33 


, -+999 
!-°-1.4 


0.53 


18 


32 


GRB050922 


1 


28.2+. 23 „ 2 
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Exponential Cutoff Model Parameters for Flares 
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Fig. 4. — Properties of exponentially cutoff power law model spectral fits to flare data for all 
Gold flares. The index number of the flares shown on the x-axis simply refers to the index 
number for each flare shown in column 1 of Table [5j Many flares did not provide enough 
data in this energy band to lead to convergence for the cutoff energy, which is clear from the 
fact that panel 3 has many data points not shown off the top of the plot (these were set to 
the 500 keV fitting limit and their lower error bars extend into the plot). 
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Table 5. Properties of exponentially cutoff power law spectral fits to Gold flares 



Index Number 


GRB 


Flare 


N H {W" J cm ') 


Photon Index 


E cut (keV) 


2 


DOF 


1 


GRB050219 


1 


40.2+ 1"' 6 , 
— 10.2 


2.68+0" 
— 0.68 


999 


1.12 


37 


2 


GRB050502 


1 




1.97+0-i 4 


7-2+ 4 n 
— 2.0 


1.36 


327 


3 


GRB050502 


3 


10.8+fl 


2.11+ -' 24 


999 


0.86 


30 


4 


GRB050607 


2 


22-B + g4 


2 - 40+ oS 


999 


0.87 


32 


5 


GRB050712 


1 


23.9+!i 


2.10+0-18 


999 


1.60 


56 


6 


GRB050712 


2 


20.2+3-0 


3.06+0-^ 


999 


1.09 


17 


7 


GRB050713 


1 


46.2+^ 




8 - 8+ -6 9 l 


1.02 


187 


8 
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2 




o oq + 0.65 
98 


999 


1.21 


47 


9 
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1 


23.6+ 6 ,f- 7 . 


i.i8+;-^ 
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55 


10 
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0.32 


54 


11 
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1 
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1.37+0-1 4 


8.7+ 4 1 
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12 
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0.95 
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13 
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21 


14 
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2 


1.6 + i 9 6 7 


0.63+. 2 - 25 , 
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0.90 


36 


15 
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1 


10.7+j!l 
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32.2+ 32 - 2 , 


1.08 


57 


Ki 
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2 
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— 0.5 


1.66+0'! 

— 0.05 


999 


0.94 


186 


17 
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3 


5.8+™ 


1 9 i+0.08 


999 


0.94 


105 


18 
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. ,„ + 0.'62 
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— 2 1 


0.98 


80 
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0.93 
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20 
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5 
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2.23+0-' 24 


999 


1.45 


33 
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999 
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Power Law + Blackbody Model Parameters for Flares 
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Fig. 5. — Properties of blackbody plus power law model spectral fits to flare data for all 
Gold flares. The index number of the flares shown on the x-axis simply refers to the index 
number for each flare shown in column 1 of Table EJ 
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Table 6. Properties of blackbody plus power law spectral fits to Gold flares 



Index Number 


GRB 


Flare 


JV H (10 20 cm- 2 ) 


Photon Index 


kT (keV) 




DOF 


1 


GRB050219 


1 


56.8+29 1 


2.34+°" 
— 0.81 


0.1+°,! 


0.96 


36 


2 


GRB050502 


1 


6.4+i'i 


2.06+.°;°! 


0.3+° : S 


1.29 


326 


3 


GRB050502 


3 




1 .86~^"q 


0.2+°| 


0.84 


29 


4 


GRB050607 


2 


26.0+ 2 . 5 s 1 




l.l+° 9 


0.89 


31 


5 


GRB050712 


1 


31.2+ 27 : 6 


2.52+ 184 


32.4+lf 3 6 


1.52 


55 


6 


GRB050712 


2 


19.5+1° 7 3 


3.04+° : | 


0.0+ 20 ° 


1.16 


16 


7 


GRB050713 


1 


41.6+^4 


, q , + 0.34 




1.01 


186 


8 


GRB050713 


2 


117.6+^1 


3.71+^ 


0.1+°° 


1.03 


46 


9 


GRB050716 


1 


31 g + 255.4 
— 31.3 


1.56+!' 36 , 
— 1.5b 


7.2+ 7 | 


0.29 


54 


10 


GRB050716 


2 




3 29 + 2-98 

— U. 78 


o.3+„ : ! 

— 0.1 


0.34 


54 


11 


GRB050724 


1 


41.7+|| 


1.58+°™ 


0.7+°-; 


1.01 


328 


12 


GRB050724 


2 


46.2+^1 


2.69+oS 


0.4+°' 3 


0.96 


52 


13 


GRB050724 


3 




i.3i +< ;- 2 = 


0.7+ -l 


1.29 


20 


14 


GRB050726 


2 


2-3±2 7 3 5 


2.03+ 3 ' 1 1 | 


0.4+° 3 


0.90 


35 


15 


GRB050730 


1 


3.2+| 2 


1.29+°^ 


0.4+2'j 


1.00 


57 


16 


GRB050730 


2 


12.6+ 2 ' 7 




2.8+ 22 8 1 


0.90 


185 


17 


GRB050730 


3 


9.6+|"° 




2 - 1+ oe 


0.91 


104 


18 


GRB050730 


4 




1.75+°' 32 
— 0.35 


o.4+° : ; 

— 0.1 


0.95 


79 


19 


GRB050802 


1 


6.2+ 29 ' 3 


2 - 45 + 0?8 


0.7+°-| 


1.00 


28 


20 


GRB050803 


5 


40. 0+ 32 -| 


"• 71 + u : 69 


0.3+°,'! 


1.16 


32 


21 


GRB050803 


6 


198.7+^1 


S.44+ 229 


0.1+°° 


0.98 


16 


22 


GRB050820 


1 


11.1+ 2 -! 


0.84+°°* 


o.o+o° 


1.12 


200 


23 


GRB050822 


1 


2.3+ 2 f 


1.32+°|° 


0.4+"- 2 

— 0.1 


0.39 


26 


24 


GRB050822 


2 


19.2+li? 


2 . 75 + 0.42 


o.i+°! 


1.09 


29 


25 


GRB050822 


3 


34.3+f ° 


4.57+j'ff 


o.o+g-g 


0.94 


16 


26 


GRB050904 


1 


14.8+ 2 '| 


1.79+ 09 


°-° + oo 


0.99 


180 


27 


GRB050904 


4 


10.4+^'° 


2.29+°'| 


l-° + 02 


1.29 


36 


28 


GRB050904 


5 


5 2+ 6 ' 5 
0-2_4.2 


1 50+ 34 


3+ 01 


0.90 


24 


29 


GRB050904 


6 


6 0+ 2 6 

D.U_ 2 g 


1 97+ 27 
-0.22 


31 ()+ 310 


1.06 


20 


30 


GRB050904 


7 


7.6+ 2 - 3 


1 96+ 20 
i.ao_ 24 


9+°' 9 
u - -0.9 


0.91 


22 


31 


GRB050916 


1 


84.1+SS;! 


, 00 + 137 
1JO _2 .69 


5 +199.4 
u '°-0.5 


0.49 


18 


32 


GRB050922 


1 




1 SR + 1 ' 93 
1 ' 86 -2.11 


2+ 01 


0.99 


97 


33 


GRB050922 


2 


20.0l 22 j 


2 65+ 127 


199 3+°' 7 


0.97 


43 


34 


GRB050922 


3 




2 3Q+0.24 


5+°- 2 


0.77 


114 


35 


GRB051117 


1 


14.3+^ 


, 77 +0.11 
1 -' '-0.12 


4+°- 3 
u *-0.1 


1.11 


340 


36 


GRB051117 


2 




2 06+ 09 


3+°'° 


1.01 


316 


37 


GRB051117 


3 


xl - D -3.0 


2 14 +0.16 

^■^-o.is 


3+°- 2 


0.99 


179 


38 


GRB051117 


4 




1 99+ 33 
i,yt, -0 .12 


„ , + 199.7 
ucl -0.3 


1.04 


224 


39 


GRB051117 




8.4t 2 ;? 


2 04+ 17 


2+°-° 


1.08 


182 


40 


GRB051117 


6 


15 l+ 2 '° 


2 20+° 12 
^ u -0.13 


5+°' 3 


1.03 


263 


41 


GRB051117 


7 


10.6i 2 '| 


2 16 + 012 


3+ 01 


1.11 


221 


42 


GRB051227 


1 




1 24+ 94 
^ -0.67 


7+°- 7 
u - ' -0.2 


0.95 


22 


43 


GRB060111 


1 


„. 7 +26.1 
i,0 -'-11.5 


9 77 +0.56 
■^■"-0.50 


3 +1 -° 


0.99 


116 


44 


GRB060111 


2 


45.0±i|;l 


, Sfi + 30 
2.Sb_ 28 


1+°° 


0.87 


74 


45 


GRB060111 


3 


30.5+j;; 7 


9 ,0 + 0.09 
2iS -0.09 


1+°° 


0.93 


295 


46 


GRB060124 


1 


i9-4i : 4 5 


, 22 + 01 
1 " ! ' ! -0.01 


o o+°-° 
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0.95 
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47 


GRB060124 


2 
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1.10 
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is skewed to the positive values of Ax 2 . One way to quantify this is to compare the number 
of GRB flares with a large A% 2 for both the real data and the simulated power law data. 
For the simulated power law data, there are 5 events from 1000 with a A% 2 > 9.0, so one 
would expect 0.23 flares in the sample of 47 Gold flares to have a A% 2 > 9.0 by chance, if 
the real data had been drawn from a power law distribution. For the real data, there are 9 
GRB flares from the sample of 47 Gold flares that have a A% 2 > 9.0. Therefore, it is unlikely 
that this sample was drawn from a simple power law spectral distribution, and it is clear 
that Band functions sometimes provide a superior fit. However, it is worth mentioning that 
a power law can provide a reasonable fit to many flares. 



4.3. Fluence of Flares 

The fluence of a flare is defined as the flux of the flare, found using the spectral fits 
described above, integrated over the duration of the flare from t start to t stop in the 0.2-10 keV 
energy band. Table [7J shows the fluences for both the power law and the Band function fits 
to the flares. Since the spectral fits provided no compelling evidence for using the thermal 
model or the exponential cutoff model, for the remainder of this paper, we restrict ourselves 
to the standard spectral models used for GRB afterglows and prompt emission, namely the 
simple power law and the Band function. The quoted error bars are la, and they include 
the error due to the uncertainty in the underlying lightcurve contribution to the fluence. 
In some cases this latter source of error is large and dominates the error from the spectral 
fit itself. The calculated fluence values for all flares (not just the Gold flares) have been 
reported, even in cases for which there are very few degrees of freedom. As a result, some of 
the fluence values are poorly constrained, as reflected by the error bars. 

These reported fluence values do not include the contribution from the power law compo- 
nent of the spectral model that was used to approximate the underlying afterglow lightcurve 
contribution to the flare spectrum. In other words, the fluence values in Column 6 of Tabled 
include only the contribution from the Band function component of the spectral model fit to 
the flare data, excluding the frozen power law component from the underlying component. 
The fluence values in Column 3 of Table [7J include only the contribution from the unfrozen 
flare power law component of the spectral model fit to the flare data, once again excluding 
the underlying component. This is an important point to stress since most previous papers 
that quote a fluence for flares actually quote the entire fluence under the lightcurve. This 
practice is misleading because the underlying afterglow lightcurve sometimes contributes a 
large, and difficult to constrain, fraction of the total fluence. The reported fluence values 
are also corrected for effects due to incomplete light curves for some flares. For instance, if 
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the tail end of a flare happened to be interrupted by an orbital gap or the South Atlantic 
Anomaly, the incomplete flare light curve would be extrapolated until it intersected with the 
extrapolation of the underlying light curve. In these cases, the error on this extrapolation 
was factored into the error on the fluence. 

The fluence calculations were done in the observed XRT energy band, which is 0.2-10 
keV. Since the X-ray flares emit the bulk of their energy in this band, we consider this to 
be the most reasonable approach. This fact is supported by the Band function spectral fits 
that result in X-ray peak energies, and it is supported by the fact that the X-ray flares 
are typically weak to undetectable by higher energy and lower energy instruments such as 
BAT and UVOT. We computed the effect of extrapolating the typical flare spectrum into the 
higher energy band typically reported for Swift-BAT bursts. The median spectral parameters 
for Gold flares that had a reasonable spectral fit were a = 1.06, j3 = 2.21, and E = 1.02 
keV; where these refer to the Band function lower energy photon index, higher energy photon 
index, and e-folding energy, respectively. From these values, it was found that extending the 
energy range from 0.2-150 keV added only 1.4% to the fluence relative to the reported 0.2-10 
keV value. This is insignificant compared to the error bars. 

The overall distribution of flare fluences is shown in Figure [7] for the absorbed power 
law model and the absorbed Band function model. Once again, these fluences are for just 
the flare component. The two distributions plotted on the right are for the 47 Gold flares 
that have > 15 DOF in the spectral fit and xled < while the two distributions plotted 
on the left include all spectral fits. Fluence derived from both power law fits (top) and 
Band function fits (bottom) are shown. The mean 0.2-10.0 keV fluence (unabsorbed) of 
our sample of flares, derived using Band function fits, is 2.4 x 10~ 7 erg cm -2 . There is no 
evidence of a bimodal distribution, which might arise if flares came from multiple processes. 
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Table 7. Fluence of Flares 



GRB 


Flare 


7 2 
.FZuence(10 erg cm ) 


2 

"red 


DOF 


7 2 
FZuence(10 erg era ) 


2 

*red 


DOF 






Power Law 




Band Function 






GRB050219 


1 


70 + 999-00 


1.09 


38 




1.12 


36 


GRB050406 


1 


n 21 +°- 02 

u ' !1 _0.07 


1.78 


7 




2.43 


6 


GRB050421 


1 


n 91 +0.17 
u - z -999.00 


0.53 


10 


0.28+"; 27 


0.64 


8 


GRB050502 


1 




1.41 


328 


8 -3°±0.44 


1.20 


326 


GRB050502 


2 


24 +017 
u - z -0.07 


1.21 


8 


o-ieto.le 


1.54 


6 


GRB050502 


3 


u - y -999.00 


0.83 


31 


°-8itS:3 


0.85 


29 


GRB050607 


1 


n on + 999.00 
u - z -999.00 


0.73 


9 


0.20+."; 44 


0.83 


7 


GRB050607 


2 


1 no + 999.00 
j..uy_ 999 |M| 


0.84 


33 




0.87 


31 


GRB050712 


1 


, e i +999 .00 
x>t>x -9»0.llll 


1.57 


57 




1.47 


55 


GRB050712 


2 


n 4Q + 999.00 
— 999.00 


1.03 


18 


2fi +0 ' 99 
u -- !D -2.44 


0.98 


16 


GRB050712 


3 


oe+999.00 
u -°°-999.00 


0.73 


10 


1 « + 0- 18 
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Fig. 6. — Histogram of A% 2 between the power law fits and Band function fits for all 
Gold flares. The histogram represents the real data, while the overlayed line represents the 
distribution of simulated power law spectra subjected to the same fitting procedure. 
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5. Flare Fluence versus Prompt Fluence 

For the purpose of comparison, the distribution of prompt emission fluences for this 
sample of 33 GRBs is shown in Figure [8) The mean prompt fluence is 2.4 x 10~ 6 erg cm -2 , 
with a standard deviation of 2.5 x 10~ 6 erg cm" 2 . The flare fluence in the 0.2-10 keV 
band, where its energy peaks, is approximately a factor of 10 less than the mean fluence of 
the prompt GRB emission (calculated in the 15-150 keV band), measured by Swift-BAT. 
However, the distributions do overlap. In at least one case (GRB 050502b), the fluence in 
a sin gle flare matches the fluence in the prompt emission from the GRB that spawned that 



flare (IFalcone et al. II2006I ). The flare fluence in the X-ray band (0.2-10 keV) is plotted as a 
function of the prompt emission fluence in the 15-150 keV band in Figure The instrument 
bands and the peak energies of the flares and the GRB prompt emission, respectively, have 
determined the bands over which we have evaluated this fluence. This is the most reasonable 
approach in the absence of more refined measurements of all spectral parameters. However, 
it should be noted that a correction to bolometric fluence could add significant fluence from 
lower energies (below 0.2 keV for flares and below 15 keV for GRB prompt emission). Based 
on the reasoning in the previous section, it is clear that extending the energy band to higher 
energies produces only insignificant effects. 



6. Flare Properties versus Underlying Afterglow Properties 

Based on rapid temporal properties, repetitive flares, and spectral changes during flares, 
past studies of individual flaring GRBs have argued that at least the flares in question were 
due to internal GRB engine properties, as opposed to afterglow related processes. This idea 
is further strengthened by the rapid rises and decays seen in the sample of flares presented 
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Fig. 7. — Unabsorbed 0.2-10.0 keV fluence distribution of flares. The two panels on the left 
are for all flares that had a convergent spectral fit. The two panels on the right are for Gold 
flares that have > 15 DOF in the spectral fit and Xred < 1-5- Fluence derived from both 
power law fits (top) and Band function fits (bottom) are shown. 
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Fig. 8. — Prompt emission 15-150 keV fluence distribution of GRBs that are in this sample 
of flaring GRBs. 
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Fig. 9. — Flare fluence in the 0.2-10 keV band (derived using a Band function) plotted as a 
function of the prompt GRB fluence in the 15-150 keV band. 
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in Paper I. In Figure [101 we compare the photon index of the power law fit to the underlying 
afterglow data to the power law photon index of the flare. The flare power laws have a wider 
distribution of spectral indices than that of the underlying afterglows. 



7. Temporal Evolution of Flare Properties 



From studies of a few individual GRBs with several strong flares, temporal evolution 



of spectral properties has 
been see n in GRB 050406 ( 



jeen explored. Spec tral evolution in individual brig 



Romano et al. Il200ffl. GRB 050502B (IFalcone et al. 



050607 flPagani et al. 1120061 ). GRB 060714 (Morris et al. 1120071 ). GRB 050822 flGodet et al. 



at flares has 



200 61), GRB 



20071 ). as well as several others. In all cases, spectral hardening was observe d at the onset 



of the flare, followed by spectral softening as the flare peaked and decayed. iKrimm et al. 



( 120071 ) discuss temporal evolution of spectral properties from flare to flare in a sequence of 
flares seen in GRB 060714, whic h show a decrease in E pea k as a function of the time of the 
flare. iButler fc Kocevski I (120071 ) has observed the same trend in his study of several bright 
X-ray flares. 

In Figure [TIJ we investigate this E pea k versus time relationship for this large sample of 
flares. We have plotted the values for all of the Gold flares that have a known redshift. The 
time axis is the flare time relative to the GRB prompt T in the burst reference frame. The 
results have not been scaled by the prompt GRB E pea k or relative to one another in any way. 
As a result, this is merely a test of whether or not an absolute relationship exists, independent 
of the individual GRB or flare parameters. There is no clear overall relationship present in 
the data. Due to the potential for unknown scaling from burst to burst, the spectral softening 
apparent in individual bursts would not necessarily be apparent when an ensemble of flares 
from many GRBs is plotted together, as we have done in this case. 



In Figure [12], we attempt to investigate the dependence of total energy release on flare 
time for this sample of flares. Once again, we have restricted ourselves to plotting the values 
for all of the Gold flares that have a known redshift. The time axis is the flare time relative 
to the the GRB prompt To in the burst reference frame. There is no clear relationship in 
the data when viewing all of the flares as a single sample as we have done here. Of course, it 
is still possible that such a temporal relationship exists for the flares in an individual GRB, 
but a scaling factor (likely to be dependent on prompt GRB parameters) would need to be 
applied to each GRB and the associated set of flares to see this effect when plotting flare 
parameters from many GRBs. 
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Fig. 10. — Comparison of the power law spectral photon indices for the underlying afterglows 
and for the flares 
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Fig. 11. — Redshift corrected peak energy of Gold flares as a function of rest frame flare time 
relative to prompt T . This plot contains all flares irrespective of (and unsealed for) prompt 
emission E peak . 
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8. E 



peak 



versus E, 



Although it is clear that the flares typically peak in the X-ray band, the limited band- 
width of the study presented here (0.3 to 10 keV) makes it difficult to constrain their peak 
energies well. However, this is a topic of considerable importance since prompt GRB emis- 
sion has shown evidence of an empirical relationship between the peak energy of the spectral 



: emission ( 


Ghirlanda et al. 


2005 


Thompson et al. 


2007 


). So, in 



have attempted to explore this relationship by looking at the relationship between the Band 
function E peak and E iso in this band. E peak is the peak energy of the Band function spec- 
trum for the flare, which corresponds to (2 + a)Eo, where Eq is the e-folding energy that 
is obtained from a Band function spectral fit. In this work, Ej SO is defined as the isotropic 
equivalent energy released during the GRB flare in the 0.2 keV to 10 MeV band, assuming a 
Band function spectrum derived in the 0.2 to 10 keV band, where the observed flare spectra 
typically peak. This is calculated as: 



E, 



a+z) 



obs\ 



(5) 



where S t> s is the unabsorbed fluence seen by the observer in the 0.2-10 keV band, z is the 
redshift, k is the correction factor from the observed 0.2-10 keV band to the co-moving 0.2 
keV to 10 MeV band, and di um is the luminosity distance calculated using a flat A dominated 
universe with Qm = 0.3 1, fl\ = 0.69, and H n = 70 km s" 1 Mpc" 1 . The k-correction factor 



is described in detail in iBloom et al. I (120011 ) . 



In Figure [131, E iso (obtained as above using Band function spectral fits) is shown as a 
function of the redshift corrected E pea k for only Gold flares that have a redshift measurement. 
Due to the paucity of jet break measurements, we can not calculate E 7 , which corrects E iso by 
accounting for the jet opening angle, th us we can not expl o re the tighter E peafc -E 7 relationship 
reported for GRB prompt emission by iGhirlanda et al. I (120051 ) . 



While it does seem clear that the flares involve a peak energy that is significantly lower 
than the more typical hundreds of keV observed for the initial GRB prompt emission, it 
is not clear if there is a strong relationship of this peak energy with E iso , due to the large 
error bars and the limited sample. The relatively low E peak in the X-ray band is, of course, 
expected since the flares are observed as increases in the X-ray band, which are often not 
accompanied by measurable increases in other bands. Unfortunately, due to the size of the 
error bars on E peak , it is not at all clear whether there is a relationship for the flares that 
is similar to the E peak -E iso correlation found for the prompt GRB emission. The intriguing 
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Fig. 12. — E iso for Gold flares as a function of rest frame flare time relative to prompt T . 
E iso is k-corrected and is calculated in the 0.2 keV to 10 MeV band. This plot contains all 
flares irrespective of (and unsealed for) prompt emission properties. 
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Fig. 13. — Exploration of Band function fit for E peak relationship with E iso for flare emission. 
Ei SO has been k corrected into the co-moving 0.2 keV to 10.0 MeV band. Only the fluence 
from the flare itself (i.e. underlying afterglow emission subtracted) was included in the 
calculation of Ej SD . The three data points with the lowest Ej SO , plotted as x symbols, are 
from flares associated with a short burst. 
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hint of a relationship evident in Figure [13] will be explored in the future using more flares 
and broader spectral coverage. 



9. Redshift Distribution 

This sample contains 14 GRBs that have a measured redshift. The redshift distribution 
is shown in Figure [TH The mean redshift for these 14 GRBs is z=2. 6. This is consistent 



with the mean redshift of all Swift GRBs, which is between 2.5 and 2.8 (IBurrows et al. 112006 



Jakobssen et al. 1120061 1 . This shows that the flaring GRBs are not drawn from a significantly 



different redshift distribution than the overall sample of GRBs. 




Redshift 



Fig. 14. — Redshift distribution of GRBs with flares. 



10. Discussion & Conclusions 

Based on this sample drawn from the first 110 Swift GRBs, it is clear that significant 
X-ray flares are produced frequently and at late times. In this paper we have presented a 
detailed spectral analysis of 77 X-ray flares drawn from 33 GRBs. Some of these GRBs had 
many flares (we find 7 significant flares in two of these GRBs, but there are probably even 
more temporally unresolved flares leading to observed intra-flare variability), while many 
GRBs had only one or two flares. Each of the flares was treated as an individual event, 
and properties of the entire sample of flares were presented. Although several studies of 
individual flares have been published recently, this paper (along with Paper I) provides a 
systematic study of a large sample of flares. Furthermore, this paper is the first to correct 
for the effects of the underlying afterglow on the spectra of the flares. For the flares in which 
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the fluence of the underlying lightcurve is a significant fraction of the fluence of the flare, 
the effect of the underlying photons on the flare spectra is significant. 

Several spectral models were fit to each flare. These models included a simple absorbed 
power law, similar to the simple absorbed power law that typically fits the underlying after- 
glow lightcurve, as well as more complex models that are more akin to the spectral models 
that frequently provide a better description of GRB prompt emission, such as the GRB 
Band function (IBand et al. 1 1 1 9 9 3h . For some flares, both simple and complex models pro- 
vided a reasonable fit, while some flares had an improved fit by using a more complex model 
such as a Band function. It is unlikely that the complete distribution of spectra was drawn 
from a pure power law distribution. Spectra with curvature, such as the Band function, 
could be related to the instantaneous source spectrum, but it should be noted that it could 
also be caused by temporal evolution of the spectrum during the flare (this work averaged 
the spectrum over the entire flare) since some flares have shown spectral evolution in time 
( jGodet et al. 1 120071 ; iKrimm et al. 1 120071 ) . In any case, this result is similar to the results 
found for prompt emission from GRBs (IKaneko et al. 1 12006c IBand et al. I Il993l ) , in which 
power laws sometimes provided a reasonable fit to prompt emission while Band functions 
provided a better fit to the overall sample. 

It was also found that the photon spectral indices of the flare spectrum did not always 
match those of the underlying spectrum and the distributions were different from one another. 
In those cases, this indicates a different population or mechanism for the production of 
the flare photons and the underlying afterglow lightcurve photons. This result provides 
further evidence that flares result from some form of internal engine activity, particularly 
when evaluated in conjunction with the small At ft values and furt her temporal analysis 
reported in Paper I, along with previous st udies of individual bursts (IBurrows et al. 1 120071 ; 
Chincarini et al. 1120061 ; iFalcone et al. 1120061 ) . Within the context of the standard model, this 
would most likely involve very late internal shocks. These internal shocks could, in principle, 
arise from a distribution of Lorentz factors for the shells from earlier internal engine activity, 
but this seems unlikely due to the ineffic iency of the kinetic energy conversion that would 
result from these weak internal shocks (jZhang I l2006bl ; lLazzati fc Pernal 120071 ). The late 
internal shocks are probably the result of late internal engine activity. 

This sample contained 14 GRBs with a measured redshift, and the average redshift 
did not differ from the average redshift for all Swift GRBs, including those without flares. 
This implies that late flares cannot be explained merely as redshifted multi-peaked emission 
from the initial prompt GRB. This result is al so supported by individu al burst analyses 
that are corrected for redshift, such as that of ICusumano et al. I (120061 ) in which a high 
redshift burst has very late flares in the burst rest frame. All of this implies that the large 
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fluence values for flares reported in this paper (sometimes comparable to the prompt GRB 
emission, and typically ~10x less than prompt GRB) must be produced at very late times 
with peak energies in the X-ray band. GRB progenitor models must be capable of producing 
this emission within a comparable energy budget to that which was previously applied to 
only the prompt GRB emission. The fallback o f material onto the central black hole af ter a 
stellar collapse could last for long time periods (jWoosley Ill993l ; iMacFadyen et al 1120011 ) and 
lead to late internal engine activity, but the reduced luminosity of this model at late times 



means that it can not explain all flares. Sever a 



mod e ls for continued ac t ivity of the centra l 



Perna et al. I (120061): iFan et al. I (I2006f ): iDai et al. I (J2006J); 



King et al. ( 2005 ); Proga and Zhang ( 2005 ); Katz (jl997 )). These models, and others 



must be evaluated within the context of the energy budget and the spectral parameters 
presented here. 

This work has also attempted to explore the r elationship between E pea k and Ei SO for 
flares, in an effort to see if the Amati relationship (lAmati 1 120061 . 120021 ) is also present for 
flare events. This is necessarily restricted to flares for which the GRB redshifts are known 
and the spectra have enough counts to constrain the parameters. Unfortunately, this leads 
to only 18 flares (3 of which are from a short burst), and the parameters are not particularly 
well constrained, as shown in Figure [131 Although a relati onship may ex ist and there is an 
intriguing hint of a correlation similar to that reported by lAmati I (120061 ) for GRB prompt 
emission, it is impossible to come to any firm conclusion due to the limited sample and 
large error bars. By looking at more flares and analyzing data from more instruments over a 
wider energy band (thus improving the E peak constraint), this aspect of the flare study will 
be revisited in the future. 
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